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Grover's database search algorithm is the optimal algorithm for finding a desired object from an 
unsorted collection of items. Although it was discovered in the context of quantum computation, 
it is simple and versatile enough to be implemented using any physical system that allows super- 
position of states, and several proposals have been made in the literature. I study a mechanical 
realisation of the algorithm using coupled simple harmonic oscillators, and construct its physical 
model for the simplest case of four identical oscillators. The identification oracle is implemented as 
an elastic refiection of the desired oscillator, and the overrelaxation operation is realised as evolution 
of the system by half an oscillation period. I derive the equations of motion, and solve them both 
analytically and by computer simulation. I extend the ideal case analysis and explore the sensitivity 
of the algorithm to changes in the initial conditions, masses of springs and damping. The amplitude 
amplification provided by the algorithm enhances the energy of the desired oscillator, while running 
the algorithm backwards spreads out the energy of the perturbed oscillator among its partners. 
The former (efficient focusing of energy into a specific oscillator) can have interesting applications 
in processes that need crossing of an energy threshold for completion, and can be useful in nan- 
otechnological devices and catalysis. The latter (efficient redistribution of energy) can be useful in 
processes requiring rapid dissipation of energy, such as shock-absorbers and vibrational shielding. I 
present some tentative proposals. 



I. GROVER'S ALGORITHM 

Database search is an elementary computational task 
with wide-ranging apphcations. Its efficiency is measured 
in terms of the number of queries one has to make to 
the database in order to find the desired item. In the 
conventional formulation of the problem, the query is a 
binary oracle (i.e. a Yes/No question). For an unsorted 
database of A'' items, starting from an unbiased state and 
using classical Boolean logic, one requires on the average 
Q = 0[N) queries to locate the desired item. 

Lov Grover discovered a search algorithm that, using 
superposition of states, reduces the number of required 
queries to Q = 0{^/N) pj. The algorithm was originally 
proposed in the context of quantum computation, but its 
applicability has since been widely expanded by realising 
that the algorithm is an amplitude amplification process 
that can be executed by a coupled set of wave modes. It 
has also been proved that the algorithm is optimal for 
unsorted database search 

Grover's algorithm starts with a superposition state, 
where each item has an equal probability to get picked, 
and evolves it to a target state where only the desired 
item can get picked. Following Dirac's notation, the 
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starting and target state satisfy (index i labels the items) , 

W\s)\''^l/N, \m\^^5u. (1) 

The algorithm evolves |s) towards \t), by discrete rota- 
tions in the two-dimensional space formed by \s) and \t). 
The rotations are performed as an alternating sequence 
of the two reflection operators, 

Ut^l-2\t){t\ , U,^l-2\s){s\ , (2) 

[-UsUt)%) = \t) . (3) 

Ut is the binary oracle which flips the sign of the target 
state amplitude, while —Us performs the reflection-in- 
the-average operation. Solution to Eq.(3) determines the 
number of queries as 

(2g-f l)sin-i(l/\/A) =7r/2 . (4) 

(In practice, Q must be an integer, while Eq.(4) may not 
have an integer solution. In such cases, the algorithm is 
stopped when the state has evolved sufficiently close to, 
although not exactly equal to, \t). Then one finds the 
desired item with a high probability.) 

The steps of the algorithm for the simplest case, Q = 1 
and A^ = 4, are illustrated in Fig.l. 

The algorithm relies on superposition of and interfer- 
ence amongst a set of states, which are generic features 
of wave dynamics. It can be executed by any system of 
coupled wave modes, provided: 

(1) The superposition of modes maintains phase coher- 
ence during evolution. 

(2) The two refiection operations (phase changes of tt for 
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FIG. 1: Amplitude evolution in Grover's algorithm for the 
simplest case, A'' = 4, when the first item is desired by the 
oracle. Dashed lines indicate the average amplitudes. 
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FIG. 2: A system of four identical harmonic oscillators, cou- 
pled to a big oscillator via the center of mass mode. 



the appropriate mode) can be efRciently implemented. 
(Note that the N states can be encoded using logj N bits, 
but have to be reahsed as N distinct wave modes.) 

Otherwise, the algorithm is fairly robust, and succeeds 
even when: 

(a) The wave modes are anharmonic though symmetric. 

(b) The initial state is somewhat randomised. 

(c) The phase changes in the reflection operations are 
slightly different from tt, 

(d) The wave modes are weakly damped. 

The interpretation of amplitude amplification occur- 
ring in the algorithm depends on the physical context. In 
the quantum version, | gives the probability of a state, 
and the algorithm solves the database search problem. 
In the classical wave version, \A\'^ gives the energy of a 
mode, and the algorithm provides the fastest method for 
energy redistribution through the phenomenon of beats. 

The quantum version of the algorithm involves highly 
fragile entanglement, and hence very short coherence 
times. It also needs to be implemented at the atomic 
scale, which is not at all easy. On the other hand, the 
classical wave version uses only superposition, which is 
much more stable, and hence it is straightforward to de- 
sign demonstration models In the following, I de- 
scribe implementation of Grover's algorithm in a simple 
mechanical setting, using four harmonic oscillators cou- 
pled via the centrc-of-mass mode. 



II. HARMONIC OSCILLATOR 
IMPLEMENTATION 

A system of coupled harmonic oscillators is frequently 
studied in physics. It involves only quadratic forms, and 
can be solved exactly in both classical and quantum do- 
mains. Let the items in the database be represented by 
N identical harmonic oscillators. While they are oscil- 
lating in a specific manner, one of them is "tapped" (i.e. 



elastically reflected). The task is to identify which of the 
oscillators has been tapped, without looking at the tap- 
ping. The optimisation criterion is to design the system 
of oscillators, and their initial state, so as to make the 
identification as quickly as possible. 

Grover's algorithm requires identical coupling between 
any pair of oscillators. That is arranged by coupling all 
the N small oscillators to a big oscillator, as shown in 
Fig. 2. The big oscillator then is coupled to the centre-of- 
mass mode, and becomes an intermediary between any 
pair of oscillators with the same coupling. Indeed, the 
centre-of-mass mode plays the role of the "average state" . 
In this setting, elastic reflection of an oscillator imple- 
ments the binary oracle in velocity space (see Fig. 3), and 
time evolution of the whole system by half an oscillation 
period carries out the reflection about average operation. 



A. Dynamics 

The equations of motion are: 

MX = -KX + k ^{xi - X) 



mxi = —k{xi — X) 



(5) 
(6) 



The non-target (xi^t) oscillators influence the dynamics 
of the target oscillator {xt) only through the centre-of- 
mass position X = J2i=i ^il^ ■ They make up TV — 2 lin- 
early independent modes, of the form Xj^t — Xk^t, which 
decouple from xt and x. Effective dynamics is thus in the 
3- variable space {X, x,xt}, with the equations of motion: 
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FIG. 3; The binary oracle flips the sign of the target oscillator 
velocity, when its displacement is zero, (a) Just before the 
oracle: Xt = 0, Vt > 0, (b) Elastic reflection: Xt = 0, — *■ 
—vt, (c) Just after the oracle: xt = 0, vt < 0. 



The last equat ion is easily solved, with the angular fre- 
quency ujt = \Jklm,. The first two equations are coupled, 
and their eigenmodes are of the form X + Aa;. We can 
find them by requiring that 

X + Xx = Xi — hA— -A— ) 



-u;^{X + \x) 



Therefore, 
A = 
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(11) 



Then the time period T = Att/loi for the whole system. 
Time evolution for half the period reverses a;, while leav- 
ing Xt — x unchanged, i.e. it implements the operator Us 
in the velocity space. Thus Grover's algorithm is realised 
by "tapping" the target oscillator at every time interval 
M = 27r/wt. 

The above resonance criterion corresponds to: 
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(17) 



In a situation where all the non-target oscillators move 
uniformly together (i.e. all Xi^t equal, all Xi^t equal), 
the displacements are 

Xt = — sin((jj+t + 0+) H — y4_ sin(cj_t + (/)_) 



-t- AtS\n{uiitt + (t>t) , 



(18) 



Xi-tt = sin(cj+t + 0+) + iy4_ sin(cj_t + (/)_) 



^At sin(wti -I- 0t) . 



(19) 



Our experimental parameters differed slightly from these 
ideal values because of various imperfections discussed 
later. 



Results 



Let the dimensionless ratios for the spring constants and 
the masses be, Rk = K/k and Rm — M/ra. Then the 
sinusoidal solutions have the angular frequency 



IN + Rk . 
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There are two solutions to these equations, coefficients 
A± and the corresponding frequencies uj±. They satisfy 



2 , 2 2 I . . ^ + Rk 

UJ^ + U}_ = IVf 1 



UJ,LU_ = LJf 
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(13) 
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The general solution to the dynamical equations is: 

X + X±x = A±sm{uj±t + (l)±) , 
xt-x = At sm{ujtt + (j>t) ■ 

B. The model 



We constructed the simplest system, with = 4. The 
binary oracle (i.e. Ut) is the elastic reflection it — > —it 
when Xt = 0. The reflection about average operation 
imposes the constraint that the whole system must evolve 
by half an oscillation period between successive oracles. 
From the many possibilities, as our design parameters, 
we selected the convenient frequency ratios 
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The uniform superposition state corresponds to the ini- 
tial conditions: 



X , X , X, 







V 



(20) 



In this case, the big oscillator returns to its initial rest 
state after every half a period, and the first binary oracle 
is applied at t = 0. The starting phases of the solution 
vanish, — (pt — 0. The amplitudes arc 
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before the binary oracle, and 

AuJt 'iuJt 

after the binary oracle. The resultant time evolution of 
the oscillators in the position and velocity spaces is illus- 
trated in Fig. 4. It is observed that Grover's algorithm 
provides position amplification of 1.87 and velocity am- 
plification of 2. 

In actual experiment, it is much easier to start with 
the initial conditions: 



X ^ A , X ^0 , Xi^^ A , ii = 



(23) 



In this case, the velocities reach their maximum value 
after a quarter period, and the first binary oracle is ap- 
plied aXt = r/4. The starting phases of the solution are, 
4'± = 0t = 7''/2. The amplitudes are 
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FIG. 4: Time evolution of the oscillators in the position (top) 
and the velocity (bottom) spaces, when the initial state is: 
X = 0, X = 0, Xi = 0, Xi = V. Red, green and blue curves 
respectively denote the uniform superposition state (without 
any oracle), the target oscillator after the oracle at t = 0, and 
the non-target oscillators after the oracle at t = 0. 



FIG. 5: Time evolution of the oscillators in the position (top) 
and the velocity (bottom) spaces, when the initial state is: 
X = A, X = Q, Xi = A, Xi — 0. Only the behaviour after 
the oracle is shown, t' = t — T/4. Red, green and blue curves 
respectively denote the uniform superposition state (without 
any oracle), the target oscillator after the oracle ai t — T/4, 
and the non-target oscillators after the oracle at t = r/4. 



before the binary oracle, and 



D. Perturbations and stability 
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after the binary oracle. The resultant time evolution (as a 
function oi t' = t — T/4) of the oscillators in the position 
and velocity spaces is illustrated in Fig. 5. It is found 
that Grover's algorithm provides position amplification 
of 1.63 and velocity amplification of 2. 

In both these cases, maximal velocity amplification is 
achieved. On the other hand, position amplification is 
substantial but not maximal. The reason is that the bi- 
nary oracle is implemented in the velocity space, and it 
is physically not possible to implement it in the position 
space. 



Gravity: A uniform gravitational field shifts the equi- 
librium positions of the vertically hanging springs, but 
apart from that it has no effect on their dynamics. This 
is obvious from the expression for the potential energy, 
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Imperfect synchronization of the initial state: When 
the initial velocities are arbitrary instead of uniform, the 
energy amplification of the target oscillator provided by 
the algorithm is not four-fold. It is instead limited to the 
initial energy present in the {X,x,xt} modes. 
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which is still substantial for the generic situation where 
the initial it and x do not differ by a large amount. 

Imprecise reflection operations: In practice, the re- 
flection operations may not exactly implement phase 
changes of tt. Also, the measurement operation termi- 
nating the algorithm may not take place at the precise 
instant of maximum amplification. The energy amplifica- 
tion depends only quadratically on such phase errors from 
the ideal values, e.g. if the reflection phase change is n+S, 
then the loss in energy amplification is 0(sin^(5/2)). 

Inelastic reflections: When the binary oracle produces 
an inelastic reflection of the target oscillator, with a coef- 
ficient of restitution 1 — e (i.e. vt (e— l)vt), the energy 
amplification decreases by the factor (1 — e/4)^. Even in 
the extreme case of e = 1, the energy amplification is a 
sizeable (9/4)-fold. 

Spring masses: Real springs are not massless. To a 
good approximation, the energy taken up by the springs 
can be estimated by assuming a constant velocity gra- 
dient along the springs, and then absorbing the spring 
energy by altering the masses of the objects attached to 
the springs. The dominant correction is to add one-third 
of the spring mass to the objects at either end. The re- 
mainder, in our set up and on the average, amounts to 
adding one- twelfth of the mass of the small springs to the 
big mass. This prescription allows tuning of the masses of 
the oscillators after measuring the masses of the springs, 
in order to maximise the energy amplification. 

Damping: For a weakly damped oscillator (damp- 
ing force —2"fmv), its amplitude changes linearly with 
the damping coefficient 7, while its frequency changes 
quadratically [1]. 



X + 2'-fx + luqX = Q 
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(28) 



Thus small external disturbances reduce the energy am- 
plification, by a factor exp(— 47r7/a;), but have little ef- 
fect on the all important phase coherence amongst the 
oscillators that governs interference of the modes. 

Overall, we find that most deviations of the param- 
eters from their ideal values affect the performance of 
the algorithm only quadratically, and can be easily taken 
care of. Damping provides the only linear perturbation, 
which should be controlled to the best possible extent. 



III. POSSIBLE APPLICATIONS 

A variety of coupled vibrational systems with small 
damping can be made easily. They can provide either fast 
focusing or fast dispersal of energy, and can therefore be 
an important component in processes sensitive to energy 
availability. 



FIG. 6: A comb-shaped cantilever beam can be used as a 
selective switch. 



A. Focusing of energy 

Efficient concentration of the total energy of a coupled 
oscillator system into a specific oscillator can be used 
as a trigger or a sensor, where an external disturbance 
becomes the cause for the reflection oracle. The focusing 
of energy could also be useful in nanomechanical systems 
[1] where the component concerned cannot be directly 
controlled, and a possibility of using coupled cantilever 
beams as a switch is pointed out in Fig. 6. 

There exist many processes that need crossing of an 
energy threshold for completion. Their rates are typically 
governed by the Boltzmann factor for the energy barrier, 
exp(— i?barrici/fc2^)- Energy amplification can speed up 
the rates of such processes by large factors, leading to 
catalysis [4j]. 



B. Dispersal of energy 

Grover's algorithm is fully reversible. The reflection 
operators Us and Ut are inverses of themselves. So the 
algorithm can be run backwards as {—UtUs)'^\t) = \s). 
That disperses large initial energy in the target oscilla- 
tor to a uniform distribution among its partners. In the 
coupled oscillator model, the initial condition would be 
if = V and Xi^t = 0. After waiting for t — 2n/Lu, and 
then reversing it produces ii = V/2. 

This behaviour can be useful in quickly reduc- 
ing localised perturbations by redistributing its energy 
throughout the system. Instead of damping a single per- 
turbed oscillator, it is much more efficient to disperse 
the energy into several oscillators while damping all of 
them together. To illustrate the concept, consider the 
situation where all the oscillators have the same damp- 
ing coefficient 7. The normal modes in the {X,x^xt} 
space then separate the same way as in Eq.(lO), and the 
relations in Eqs.(12,13) are retained. Damping shifts the 

and 



oscillation frequencies according to uj' = y'u 
the general solution becomes: 



X + X±x = A±e-'^* sm{Lj'^t + , (29) 
xt~x ^ Ate-'^*sin(wji + 0t) . (30) 
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FIG. 7: A hierarchical system of coupled oscillators can act 
as a shock absorber. The initial impulse is assumed to be a 
local disturbance, which subsequently spreads out. 



The coupled oscillator dynamics of Grover's algorithm is 
maintained by keeping the frequency ratios unchanged, 
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2 ' ' " 2 
If the initial conditions arc chosen as 

X^O , X = , = 0, xt^V , x,^t = , (32) 

then after half an oscillation period, T'/2 = 2ti/ijj[^ 
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This results show that the distribution of energy among 
the coupled oscillators suppresses the energy of the target 
oscillator by an extra factor of 4, in addition to the usual 
damping factor for a stand-alone oscillator. It is indeed 
the maximum possible reduction in energy, combining 
both the mechanisms. 

The damping coefficient that maximises the rate of en- 
ergy loss of the target oscillator is given by 



d_ 



-47r7 / u. 
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7 = 0.213wt 



(34) 



With this choice, in time T' /2, the energy of the target 
oscillator is reduced to 1.6% of its initial value. Although 
the energy loss is more compared to a stand-alone oscil- 
lator, it is less localized because it is distributed among 
the target oscillator, its partners and the big oscillator. 
The corresponding spring constant and mass ratios are: 



Rk — 2.91 , Rr, 



4.68 



(35) 



A hierarchical system of coupled oscillators (see Fig. 7) 
can be even more efficient in dispersal of energy, by imple- 
menting the above mechanism simultaneously at multiple 
scales. The simplest choice would be to couple four small 
oscillators to a big one at every level, with appropriate 
mass, spring and damping parameters. 

C. Transfer of energy 

It is also possible to combine dispersal and concentra- 
tion operations to to transfer energy from one oscillator 
to another via the centre-of-mass mode. For example, in 
case of four coupled oscillators, initial energy in oscillator 
\ti) can be transferred to oscillator [^2) by 



iUMt,)iUt,Us)\h) - |i2> 



(36) 



In this manner, a local signal received by a large detector 
can be first dispersed over the whole system and then 
extracted at a specific location. 
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